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The Influence of Horizontal Inhomogeneity 
on Radiative Characteristics of Clouds: 
An Asymptotic Case Study 


Alexander A. Kokhanovsky 


Abstract—The paper is devoted to analytical studies of the influ- 
ence of horizontal inhomogeneity of clouds on their radiative prop- 
erties in the framework of the asymptotic radiative transfer theory. 
It is assumed that a cloud field is optically thick. Thus, only overcast 
cloud fields are under consideration. The study is based on the in- 
dependent column approximation, assuming Gamma distribution 
of cloud optical thickness (COT) in the cloud field under study. This 
paper confirms all essential findings of cloud optics, concerning the 
influence of horizontally inhomogeneous clouds on transmitted, re- 
flected, and absorbed solar light. For instance, we found the de- 
crease of light reflection and absorption (negative biases) and the 
increase of light transmission (positive biases) for inhomogeneous 
clouds as compared to the case of a homogeneous cloud field with 
the value of optical thickness equal to that of an average COT 
of an inhomogeneous cloud field. Analytical equations for biases 
of radiative characteristics (e.g., the reflection function, spherical 
albedo, transmittance, and absorptance) are derived. Also, for the 
first time, we established the relationships between biases of dif- 
ferent radiative characteristics. Analytical equations proposed are 
simple and highly accurate for optically thick weakly absorbing 
cloud fields. They can be incorporated in large-scale atmospheric 
models for the simplification of radiation blocks both for homoge- 
neous and inhomogeneous cloud fields in visible and near-infrared 
spectral regions. 


Index Terms—Clouds, radiative transfer, remote sensing. 


I. INTRODUCTION 


HE INFLUENCE of the horizontal inhomogeneity of 

clouds on their radiative characteristics is a major subject 
of modern cloud optics studies [3], [6], [7], [9], [10], [12], 
[25], [26], [38], [43], [44]. In particular, it was found that the 
horizontal inhomogeneity of clouds effects their abilities to 
absorb, reflect, and transmit solar light [11], [44]. Thus, cloud 
remote sensing techniques, based on the spectral reflectance 
method [1], [14], [18], [22], [32], [33], [40]-[42], should 
account for the subpixel cloud horizontal inhomogeneity. This 
is not generally the case so far. 

It is known that pixels with inhomogeneous clouds are darker 
than pixels with homogeneous cloud layers having the same av- 
erage optical thickness [6]. This leads to the underestimation 
of cloud optical thickness (COT) by modern satellite retrieval 
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techniques. There are semiempirical approaches to overcome 
this problem. They are based on the artificial increase of the 
reflection function measured to account for the horizontal inho- 
mogeneity of a cloud field under study. The correct magnitude 
of these adjustments, however, cannot be assumed a priori. So 
they lack a physical basis. This issue is discussed in detail by 
Pincus and Klein [37]. Also cloud inhomogeneity could lead to 
unphysical dependencies of COT retrieved on illumination and 
viewing geometry [25]. 

Another way to solve the problem is to use three-dimensional 
(3-D) Monte Carlo calculations (e.g., see [44]). However, they 
are time consuming for realistic clouds and can be used mostly 
for theoretical studies and not as a core of operational cloud 
satellite retrieval algorithms. Monte Carlo calculations have 
shown, however, that in some cases a high accuracy can be 
achieved if a 3-D cloud field is substituted by N noninteracting 
vertical columns or cells. A cloud field in each cell is modeled 
as a horizontally homogeneous plane-parallel layer of an 
infinite horizontal extent. The optical thickness (and possibly 
microstructure) of each cell varies, depending on its position 
in a cloud field. Such an approach is called the independent 
column approximation (ICA) or the independent pixel approxi- 
mation (IPA). The range of applicability of the ICA was studied 
by Davis et al. [8] and Schreier and Macke [44]. 

Effectively, the ICA reduces the 3-D radiative transfer 
problem to N standard radiative transfer problems for homo- 
geneous media. The number N can be large. Thus, the problem 
remains computationally very expensive. 

It can be simplified, however, if one applies approximate so- 
lutions of the radiative transfer problem for each cloud cell. It is 
done usually in the framework of the two-stream approximation 
[2]-[4], [36]. However, the accuracy of two-stream approxima- 
tions [29] is rather low as compared to exact radiative transfer 
calculations [16], [49]. In particular, for some cases, errors in- 
troduced by the approximation can be larger than differences 
of radiative fluxes for horizontally homogeneous and inhomo- 
geneous cloud fields itself. This approximation also does not 
allow to consider the bidirectional reflection function of clouds, 
which is routinely measured by various radiometers and spec- 
trometers on satellite platforms. With this in mind, we propose 
here to use the asymptotic equations [51] of the radiative transfer 
theory to solve each of N standard radiative problems, discussed 
above. Asymptotic formulas are much more accurate than those 
of the two-stream theory [16]. Typically, the error of asymp- 
totic equations is below 1% for cloud reflection functions with 
optical thicknesses 7 > 9 [17]. 
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The unability to deal with thin clouds is clearly a drawback 
of the approach presented here. However, the ICA itself is in 
trouble in this case due to the horizontal photon transport [38]. 
Moreover, a great portion of clouds has COTs larger than 5—10 
[50]. Thus, we believe that results presented here can be used in 
most of cases, occuring in cloudy atmospheres (and always for 
overcast conditions). 

The main idea behind this paper (apart from the theoretical 
model proposed) is to quantify the effects of inhomogeneity of 
clouds on their radiative characteristics in visible and near-in- 
frared spectral ranges, where solar light absorption by clouds is 
comparatively low. 

Our results can be also used in global circulation models 
(GCMs), which currently entirely neglect the cloud horizontal 
inhomogeneity [27], [45]. In particular, this drawback leads 
to inconsistency between the hydrologic and radiative transfer 
water content treatments in modern GCMs. 


II. COT DISTRIBUTION 


It is known that even overcast cloud fields cannot be charac- 
terized by a constant geometrical thickness H [3]. Instead, H 
varies along the cloud field. The same is true for the COT 7 and 
also for other cloud parameters and properties, including their 
reflection and transmission functions. The statistical theory is, 
therefore, a natural tool for cloud studies [28]. 

In this paper, we ignore the variability of cloud microphysical 
properties (e.g., the size and shape of particles) inside clouds. 
This allows us to concentrate on the influence of effects of COT 
variability on cloud reflection and transmission functions. For 
this, we need to make an assumption on the COT probability dis- 
tribution function (pdf) f(T). Generally speaking, the COT pdf 
is not a contiuous function, but rather a discrete one. However, 
for the sake of simplicity itis usually assumed to be a continuous 
function. Satellite and airborne measurements show that f(7) 
can be well represented by the lognormal, beta, or gamma distri- 
bution functions [33], [35]. The gamma distribution is the most 
convenient for theoretical studies. So, we will use the gamma 
distribution in this paper. Note, however, that there is a simi- 
larity between gamma and log-normal distributions [19]. Also, 
global optical characteristics of clouds are most sensitive to the 
average COT 


r= | idr (1) 


and the standard deviation 


ices y f me e 2) 


The particular choice of f(T), therefore, is of a secondary im- 
portance. 
Taking into account mentioned above, we assume that 


To 


f(r) = Dr" exp (-n7) (3) 
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Fig. 1. COT distribution f(T) at various values of p+. 


where 


D=—————_ 4 
a EGET) = 


and T (u + 1) is the Gamma function, defined as follows [23]: 


r(e) = | ate "da. (5) 
0 
It holds that 


~“ f(r)dr =1. (6) 


Broken cloud fields, which are out of scope of this paper, 
can be studied using other approaches (e.g., see Zuev and Titov 
[53]). 

Parameters 7) and p in (3) have a clear meaning. Namely, 
To specifies the mode optical thickness, where f(T) has a max- 
imum, and p is the half-width parameter. The distribution is nar- 
rower for larger p. 

Let us introduce the coefficient of variance 


al 


pr = (7) 


The parameter p, is equal to the ratio standard devia- 
tion/mean and has a clear statistical sense. It is often called the 
inhomogeneity parameter. The pdf (3) is given in Fig. | for 
various values of p;. We have for the distribution (3) 


F= (1 + 1) (8) 
m 


(9) 


Parameters 7 and p, can be defined for any distribution func- 
tion f(T). Therefore, they are used in the discussion, which fol- 
lows. 

Note that it follows from (8) and (9) 


w= pp? — 1, T = (1 — 02). (10) 
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Fig. 2. Values of p for various scenes, defined by Barker et al. [3]. 


It follows that u — 00, To —> T as pr — 0, as it should be. 
The value of 7 for warm clouds is usually in the range 5-50 [50]. 

The value of p, depends on many factors, including the 
averaging scale. Experimantal data for the value of p, of natural 
cloud fields are rare. Nakajima et al. [33] found that the value 
of p+ is in the range of 0.18—0.27 for cloud systems they studied 
(see Fig. 1). Their functions f(7) followed fairly well the log- 
normal distribution in contrast to highly skewed distributions 
with long tails for larger 7, reported by Gorodetsky et al. [13], 
King [17], and Heidinger and Stowe [15]. King [17] identified 
a multimodal nature of the distribution f(7) for the cloud 
system he studied with main modes, located at 7 equal to 7, 33, 
and 45. Our estimation of p,, obtained from data of Heidinger 
and Stowe [15] suggests that p, ~ 0.4 for their study. Data for 
p+, calculated from Barker et al. [3, Table II], are presented in 
Fig. 2. They were obtained analyzing 18 Landsat scenes, having 
2048 x 2048 (scenes 1-5) or 1024 x 1024 (scenes 6-18) pixels 
at 28.5- and 57-m resolution, respectively [3]. Pixel values of 7 
were inferred, using the standard cloud retrieval algorithm for 
plane-parallel slabs [33]. We see that the value of p, is in the 
range 0.2—0.8 for cloud fields studied. The approximate cloud 
field area ©) was 58 x 58 km? [3] for each scene. 

Note that some spectrometers on satellite platforms have the 
spatial resolution, which is close to (or even larger than) values 
of X} reported above. 

Thus, it is of importance to study the influence of the param- 
eter p, on the cloud radiative characteristics. It will be done in 
this paper assuming that p+ = 1/4, 1/3, 1/2, which cover most 
of cases presented in Fig. 2 (see also Fig. 1). 

It should be pointed out, however, that our study has only a 
limited value as far as a subpixel cloud inhomogeneity is con- 
cerned. Indeed, the current generation of satellite imaging sen- 
sors (e.g., see King et al. [18] and Nakajima et al. [34]) has pixel 
sizes of 1 km and even less than that, and therefore, the bias p 
is considerably reduced in this case. Also, the horizontal pixel 
inhomogeneity is of importance for small pixels, which is not 
accounted for in our model. On the other hand, some modern 
spectrometers (e.g., [5]) have a large spatial resolution (e.g., 
30 x 60 km”). Then our analysis is relevant not only in respect 
to mesoscale cloud properties studies but also for the estimation 
of influence of subpixel cloud inhomogeneity. 


III. VISIBLE RANGE 
A. General Equations 


The information on the function f(7) allows us to study the 
optical characteristics of inhomogeneous clouds in the frame- 
work of the independent column approximation, which neglects 
horizontal photon transport. The reflection R and transmission 
T functions of cloud fields depend on f(r). It is not correct 
to assume that R(T) and T(7) give measured values of reflec- 
tion and transmission. The relationship-measured signals with 
the statistical distribution f(7) is more involved. Also, the hori- 
zontal photon transport can complicate correspondent relation- 
ships. For large pixels (e.g., 30 x 60 km7), the effects of hori- 
zontal photon transport are largely canceled out, and measured 
values can be well represented by 


R= | R(r)f (r)dr (11) 


and 


CO 
T= l Reri (12) 
0 
where f(T) is the COT pdf. Equations (11) and (12) are of im- 
portance not only in subpixel inhomogeneity influence studies 
(for large pixels), but also in deriving mesoscale cloud radiative 
properties from measured COT fields f(T), which is of impor- 
tance, for example, for climate change studies. 

Other statistical characteristics of random fields R(T) and 
T(r), which are determined by the statistical distribution of the 
COT, can be considered as well, e.g., we have for the coefficient 
of variance 


(13) 


SIS 


Pr = 


oa fe — T)? f(T)dr. 


Here x = R, x = T or any other relevant cloud radiative charac- 
teristic. It is usually assumed in cloud properties retrievals that 
[14] 


where 


(14) 


ZH U(r): (15) 
This relation is correct only in the limit py — 0 (see Fig. 1). 
Thus, it is important to study biases 
By, =1- 2o (16) 
T 
for various cloud radiative characteristics x. This is also of im- 
portance for mesoscale cloud radiative characteristics studies. 
The main idea of this paper is to study dependencies B, and 
Px on T and p, for overcast clouds, based on the asymptotic ra- 
diative transfer theory valid as r > 9 [17]. We study a large 
number of cloud characteristics, including the plane albedo, re- 
flection function, cloud absorptance, and cloud transmittance 
(see Table I). Note that light field polarization characteristics 
can be considered in a similar way. 
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TABLE I 
DEFINITIONS OF VARIOUS RADIATIVE CHARACTERISTICS [24], [30]. HERE, T* = 1/27 Je T(vo, 0, ¢)dy, R* = 1/27 fe" R00, V, y)dy, AND Fo IS THE 


INCIDENT FLUX DENSITY. F= IS THE TOTAL OUTGOING FLUX (“++” FOR THE REFLECTED LIGHT FLUX AND ““—” OTHERWISE). I> 


= IS THE LIGHT FIELD INTENSITY 


(“+” FOR THE REFLECTED LIGHT AND “—” OTHERWISE). o, V, ARE THE SOLAR ZENITH ANGLE, OBSERVATION ZENITH ANGLE, AND AZIMUTH, RESPECTIVELY 


Radiative characterisic Symbol Definition 
Reflection function R (9,0, p) Ron 
Transmission function T (Bo, V, p) EAN 
T/2 
Reflectance (plane albedo) Ra (80) Roa = J R (0) sin 2009 
0 
n/2 
Diffuse transmittance Ta (vo) Hoste = f T*(%)sin 20d 
0 
n/2 
Global reflectance (spherical albedo) r J Ra (Vo) sin 2Wod¥o 
0 
mw /2 
Global transmittance t J Ta (do) sin 2od¥o 
0 
Absorptance Aa (vo) 1 — Ra (o) — Ta (Vo) 
Global absorptance a l-r-t 
TABLE II 


MODIFIED ASYMPTOTIC EQUATIONS OF THE RADIATIVE TRANSFER THEORY [20]. HERE, rə = exp(—y) Is THE SPHERICAL ALBEDO OF A SEMIINFINITE LAYER. 
Raco = exp(—yKo(V0)), Roo (0, 0, p) = R? Wo, 9, ¢) exp(—yu (Vo, V, ~)), u = Ko (9) Ko (V0) /R (00.0, p), £ = T1/3(1 — g)(1 — wo), 
y = 4/1 — wo) /3(1 — g), Ko(¥0) = (3/7)(1 + 2 cos Jo), v & 1.07 


Radiative characterisic Symbol Equation 
R v 0 p) — 
Reflection function R (0o, ù, p) o (Bo ) 
(Too — r) Ko (8) Ko (Bo) 
Transmission function T (00,8, p) tKo (0) Ko (vo) 
Reflectance (plane albedo) Ra (vo) Raco — (Too — r) Ko (vo) 
Transmittance Ta (vo) tKo (vo) 
Global reflectance (spherical albedo) r soho toy) exp yew expe) 
Global transmittance t ae 
Absorptance Aa (vo) 1 — Ra (vo) — Ta (vo) 
Global absorptance a l-r-t 
We also consider the relationships between biases of various where 
cloud radiative transfer characteristics (e.g., Br and Br, etc.). 3 
The general form of such relationships can be easily found, ma z} (19) 
taking into account that a specific directional cloud radiative B 


characteristic x of optically thick cloud fields (e.g., bidirec- 
tional reflection function) can be presented as a linear form (see 
Table II) of the angular averaged value of x, which we denote s 

z(T) =a + Bs(r) (17) 
where a and 8 depend on the viewing and illumination geom- 
etry (see the Tables I-III), but not on 7. Note that parameters a 
and 8 differ for different x (see Tables I-III). It follows from 
(16) and (17) that 


Bz = Bs (18) 


Equations (18) and (19) allow us to reduce the calculation of 
biases for angular depend characteristics x to the calculation of 
biases in the values of angular integrated characteristics s, which 
we can consider as an important result of this study. Numerous 
specific results of such relationships are given below. 

We also note that it follows due to a linearity of (17) 


By Px 
pitied eee ce (20) 
Bs Ps 
Therefore 
Px = ly|Ps (21) 
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with the same value of y as in (18). This means that it is enough 
to consider the relationship between biases. Relationships be- 
tween coefficients of variance are similar in form. The differ- 
ence arises only due to the fact that biases could be negative 
and positive. By definition, coefficients of variance are positive 
numbers [see (13)]. 

Note that y = 1 (and therefore By = B, and py = p,) at 
a = 0, which is the case, for example, for transmitted light (see 
Tables I-III). This means, in particular, that biases for transmis- 
sion function, diffuse transmittance, and global transmittance 
coincide in the case of optically thick cloud fields. We consider 
this point in some detail in Section II-B. 


B. Light Transmission 


The cloud transmission function does not depend on the az- 
imuth ¢ for optically thick clouds and is given by the following 
approximate equation in visible, where light absorption by water 
droplets can be neglected [21], [51] 


T(r, £n) = Ta(t)Ko(n) Talr) =t(7)Ko(£) D 
with a high accuracy. It follows approximately [46] 
3 
Ko(€) = 7 + 26) (23) 
and 
UT) = —— (24) 
Oe aia 


Here ¢(7) is the global transmittance (see Tables I-III). Here 
v x 1.07, € = cos(¥o), n = cos(#), J, and Jp are observa- 
tion and incidence angles, respectively, and g is the asymmetry 
parameter of a cloudy medium, which is assumed to be equal 
to 0.85 throughout this study. The value of T(r) is called the 
diffuse transmittance (see Tables I-III). Using results given in 
Section III-A, we easily obtain 


T =tKo()Ko(n) (25) 
Br = Br, = B, (26) 
PT = PT, = Pt: (27) 


We see, therefore, that the bias of the transmission function 
Br (and also the diffuse transmittance bias Br,) is equal to the 
bias for the total transmittance B+. The same is true for coeffi- 
cients of variance. This is due to the proportionality of functions 
T and Ty to the value of ¢ and the independence of the function 
Ko(€) on r. 

We have for t, B;, and p; [see (12), (13), (16), and (24)] in 
the case of the Gamma distribution (3) 


t =v! (v, v, z) (28) 
Purr Gee Be A we 
a= y OA =i (30) 
Here 
nab EREE on 


is the confluent hypegeometric function of the second kind [23] 
and 


4vu 


cr) A 


v=1+u= p7”. 

We see, therefore, that values of t, B+, and p; are determined 
by two parameters, namely v and z, which are related to p, and 
7. Using the asymptotical behavior of the function U(e, d, z) at 
small z, we obtain as T — oo 


4 
t= 33 
(1 — gy) G 
Beep; (34) 
Pr (35) 


Therefore, the problem of finding biases for the transmission 
function T and the diffused transmittance Ty is reduced in the 
case studied to the same problem but for the global transmit- 
tance t. This allows to avoid geometric considerations (e.g., the 
account for angles Yo, 9, and p). Also, it means that biases in the 
transmission J are angular independent. They are determined 
only by values 7 and p, in the Gamma-independent column ap- 
proximation studied here. 

Functions B;(7) and p;(7), obtained from (29) and (30), are 
presented in Fig. 3(a) and (b) at p- = 1/4, 1/3, 1/2. We see that 
both B;(7) and p; (F) increase with 7. They tend to asymptotic 
limits, given by (34) and (35), respectively, as 7 — oo. 

The bias B; is positive. It means that t(7) < t [see (16)]. One 
concludes from Fig. 3(a) that B; can reach 20% at 7 = 100 and 
pr = 0.5, which is not a small number as far as climate change 
and cloud remote sensing algorithms are of concern. 


C. Light Reflection 


Let us consider now quantities (11), (13), and (16) for the 
reflected light. For this we note that the reflection function of 
optically thick clouds is given by the following expression in the 
visible, where light absorption by clouds can be often neglected 
[20], [51], [52]: 


RE, N, Ps T) = R(E, n, p) E t(7)Ko(€)Ko(n). 


By definition, R3(€,7,~) = R(E,n,~,T — 00) does not 
depend on 7. So we have taking into account results given in 
Section II-A 


(36) 


R(E,n, 9,7) = R% (En, e)(1 — ut) (37) 
Br = - MB; (38) 
pr =|M]|p: (39) 
where 
ut 

ME 1 — ut (40) 

Ko(€)Ko(n) 

Sis 41 

P(E.) ve 


We see, therefore, that the problem of finding biases for the 
reflection function is reduced to the problem studied in Sec- 
tion III-B. So, having information on 7, B;, and p+, one can 
easily derive R, Bpr, and pr. For this, also, we need to know 
the function R? (£, n, p) = R(E, n, p, T —> 00) (e.g., see [31]). 
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Fig. 3. (a) Dependencies B,(T) > 0, B,(T) < 0 (in percent) at p- = 1/4, 


1/3, 1/2. (b) The dependence p;(F) at p+ = 1/4, 1/3, 1/2. 


This function practically does not depend on the cloud micro- 
physical parameters and can be considered as a universal func- 
tion in most of cases [19], [20]. 

We present the amplification coefficient M, given by (40), in 
Fig. 4. It increases from 0 to 1.5 while b = ut changes from 0 
to 0.6. For optically thick clouds, ¢ is usually a small number, 
and M < 1, which means that Br and pp are smaller than Br 
and pr, respectively. Clearly, for semiinfinite clouds, we have 
Br = Oand pr = 0. 

The function u(v%) at nadir observation [see (41)] is 
given in Fig. 5, where we used the analytical expression for 
R? (Vo, V = 0), given by Kokhanovsky [20]. The value of 
u [and therefore M; see (40)] generally decreases with the 
solar zenith angle. Thus, low sun will lead to smaller values of 
|Br| and pp [see (38) and (39)]. This can produce unphysical 
dependencies 7(W%9) in operational cloud retrieval algorithms. 

We see that M > 0. This means that biases for reflection and 
transmission have different signs. They are positive for the trans- 


Fig. 4. Dependence of the amplification coefficient M on the parameter b = 
ut [see (40)]. 


o 20 40 60 80 
solar angle, degrees 


Fig. 5. Dependence (Vo), obtained from (41). 


mission, but negative for the reflection. Inhomogeneous clouds, 
therefore, reflect smaller amount of light as compared to homo- 
geneous clouds with the value of r = 7. The opposite is true 
for the transmission, which increases for inhomogeneous cloud 
fields [52]. 

The dependencies Bp(T) and pr(T) at pr = 1/4,1/3, 1/2, 
obtained from (38) and (39), are given in Figs. 6 and 7 (Jo = 
60°, V = 0°). As expected, |B] decreases with 7 . It increases 
with p7. However, the bias is not negligable for smaller T, e.g., 
it is 8.5% at T = 10 and p; = 0.5 (see Fig. 6). Such biases 
cannot be ignored in satellite cloud retrival algorithms, based 
on the passive remote sensing techniques. 

Let us briefly consider the case of the diffuse reflectance and 
spherical albedo (see Tables I-III). Then we have 


Ra(€) =1—tKo(€) (42) 
Br, = —NB, (43) 
pr = |N| o (44) 
where 
Ko(&)t 
N = ———. 45 
1 — Ko(é)t = 


Here we accounted for the fact that Ry, and roo (see Table II) 
are both equal to one for nonabsorbing media [19]. We see again 
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TABLE III 
AVERAGE RADIATIVE CHARACTERISTICS (A goo = 1 — Racos doo = 1 — Træ) 


R= Ræ-~ (Too — F) Ko (8) Ko (vo) 
Ra= Rao- (Too — F) Ko (vo) 

T= {Ko (0) Ko (vo) 

Ta= tKo(vo) 

Aa= Ado + (@ — ao) Ko (Vo) 

a= 1-F-t 


that the multiplier N > 0 as in the case of reflection function. 
It follows for the spherical albedo 7 and values B,. and p,. (see 
Tables I-IID 


r=1-t B =-LB; pr= Lp (46) 


where 


(47) 


The dependence B,.(7) is given in Fig. 3(a). We see that it 
is negative and |B,.| < |B,|. Thus, effects of clouds inhomo- 
geneity are more important in the transmission than in the re- 
flection for a special case of overcast cloud fields considered 
here. 

We undeline once more that biases of different radiative char- 
acteristics are not independent. The analytical relationships be- 
tween them are given above. 


IV. NEAR-INFRARED RANGE 


Let us consider now the near-infrared range, where cloud par- 
ticles absorb incident solar light. The exact asymptotic theory 
leads to quite complex expressions for reflection and transmis- 
sion functions, if the value of the single scattering albedo wo 
is arbitrary [51]. However, wo is close to one for clouds in the 
near-infrared region of the electromagnetic spectrum [19]. It al- 
lows to substitute exact asymptotic formulas by more simple 
relations almost not loosing the accuracy [20], [52]. These re- 
lations will be used in this section to consider the case of ab- 
sorbing inhomogeneous clouds. They are given in Table II. Cor- 
respondent equations for various average radiative transfer char- 
acteristics are presented in Table II. 

We see that radiative transfer characteristics of inhomoge- 
neous clouds depend on Roo, Too, T, t, and Ko(Vo). Functions 
Roo; Ko(ðo), and roo are given in Tables I-III and 


r= f eer Fe f e aS 
(0) 0 


where [20] 


e sinh(x + vy) exp(—y) — sinh(y) exp(—x — y) (49) 
sinh(a + vy) 
and 
E sinh y 
iN sinh(x + vy)” ey 
Here v & 1.07 as specified above and 
Sts Vea. aig 


3(1 — g) 


It is easy to show in a full anology with results obtained above 
that 


Br=PB, Br,=QBr Br=Br,=B, (5) 


where 
Tu Q= FKo(Vo) 

1- (Too - T)Ko(ðo) ` 
Expressions similar in form hold for coefficients of variance 
[see (21)]. 

Let us introduce also the absorption function A = 1- R-T, 
the absorptance Ag = 1 — Ra — Ty, and the global absorptance 
a= 1 — r — t. The value of Ag gives a fraction of solar energy 


absorbed by a cloud field for a given sun zenith angle. Then we 
have 


(53) 


A=1-R-T Ag=1—-Ra-—Ra @=1-7F-#t (54) 
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Fig. 8. (a) Dependence B,(7) (in percent) at the probability of photon 


absorption 3 = 0.01, the asymmetry parameter g = 0.85, and p- = 1/4, 
1/3, 1/2. (b) Dependence B,(7) (in percent) for the same conditions as in (a). 


and 
Ba=VB, Ba, =WBa (55) 
where [see (19)] 
a a 
See a ae (39) 
Here 
ics l= Reo - oo Ko(€) Ko(n) 
Ko(€)Ko(n) 
1- Rä = Aco K(€) 
O= 57 
Ko(€) oe 


and aœ = 1— ræ. 

Thus, the biases Br, Br, B4, and correspondent coefficients 
of variance pr, pr, and p4 are easily obtained if the values of 
Br, Bi, Ba, Pr» Pt, and pa for absorbing inhomogeneous clouds 
are known. 

The results of the numerical calculation of values B, and Ba 
using (16) and formulas given in the Tables I-III for absorbing 
clouds are presented in Fig. 8(a) and (b). It was assumed that 
B = 1— wo = 0.01 and g = 0.85. It follows that | B,.| decreases 
with 7 and B, < 0. This corresponds to similar results given in 


Fig. 3(a) for nonabsorbing clouds. The dependence | B,(7)| has 
a maximum, which is located approximately at 7 = 30 for the 
case studied. The fact that B,(T) < 0 points out that inhomo- 
geneity of clouds leads to the decrease of the cloud absorption 
as compared to the case of homogeneous clouds with the same 
average optical thickness (@ < a(T)). 


V. CONCLUSION 


The study of the influence of the horizontal inhomogeneity 
on cloud radiative characteristics are usually performed using 
Monte Carlo techniques, which are computationally very expen- 
sive. They are related only to the specific cases computed. Other 
approaches (e.g., the two-flux approximation) are characterized 
by a poor accuracy and are limited to the calculation of light 
fluxes. 

In this paper, we propose the use of the asymptotic equations 
of the radiative transfer theory, which are valid at 7 > 9, in com- 
bination with the independent column approximation to study 
the influence of the horizontal inhomogeneity on the radiative 
properties of overcast cloudiness. 

Final formulas (e.g., see Tables I-III) are simple and allow for 
a rapid estimation of cloud inhomogeneity effects in the asymp- 
totic limit of large 7. In particular, we confirm that cloud inho- 
mogeneity leads to the increase of transmission and the decrease 
of reflection and absorption as compared to the case of homoge- 
neous clouds with the same average optical thickness. This may 
complicate so-called cloud absorption paradox [39], [47], [48] 
and have an importance for climatological studies. 

We find that biases for all radiative characteristics studied can 
be expressed via the bias B; (or B,.) in the case of nonabsorbing 
(e.g., in the visible range) optically thick clouds. For absorbing 
clouds (e.g., in the near-infrared range), we have established re- 
lationships between biases of angular integrated and directional 
quintities. 

The procedure for calculation of averages given here can be 
used to improve radiation blocks in modern global circulation 
models for the case of overcast cloud fields. 
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